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Abstract 

Assuming that the mechanism proposed by Gell-Mann and Hartle works as a mech- 
anism for decoherence and classicalization of the metric field, we formally derive the 
form of an effective theory for the gravitational field in a semiclassical regime. This 
effective theory takes the form of the usual semiclassical theory of gravity, based on 
the semiclassical Einstein equation, plus a stochastic correction which accounts for the 
back reaction of the lowest order matter stress-energy fluctuations. 



*Institut de Fisica dAltes Energies (IFAE) 



I. Introduction 



In the semiclassical theory of gravity, the gravitational field is treated classically, but the 
matter fields are quantum. The key equation of the theory is the semiclassical Einstein 
equation, a generalization of the Einstein equation where the expectation value of the stress- 
energy tensor of quantum matter fields is the source of curvature. 

One expects that semiclassical gravity could be derived from a fundamental quantum 
theory of gravity as a certain approximation, but, in the absence of such a fundamental 
theory, the scope and limits of the semiclassical theory are not very well understood. It 
seems clear, nevertheless, that it should not be valid unless gravitational fluctuations are 
negligibly small. This condition may break down when the matter stress-energy has appre- 
ciable quantum fluctuations, since one would expect that fluctuations in the stress-energy of 
matter would induce gravitational fluctuations |TJ . A number of examples have been recently 
studied, both in cosmological and flat spacetimes, where, for some states of the matter fields, 
the stress-energy tensor have significant fluctuations @. To account for such fluctuations, it 
is necessary to extend the semiclassical theory of gravity. 

To address this problem, or analogous problems in quantum mechanics or quantum field 
theory, different approaches have been adopted in the literature. The present paper attempts 
to unify, at least conceptually, two of these approaches in a formal derivation of an effective 
theory for the gravitational field in the semiclassical regime. The common feature of these 
two approaches is the idea of viewing the metric field as the system of interest and the matter 
fields as being part of its environment. This idea was first proposed by Hu [f§ in the context 
of semiclassical cosmology. Both approaches make use of the influence functional formalism, 
introduced by Feynman and Vernon |4| to deal with a system-environment interaction in a 
full quantum theory. In this formalism, the integration of the environmental variables in a 
path integral yields the influence functional, from which one can define an effective action 
for the dynamics of the system [g g 0, g, g 0, . 

The first of these two approaches has been extensively used in the literature, not only in 
the framework of semiclassical cosmology [7|, g, |T|, |TJ], |T5l [LlJ , but also in the context of 
analogous semiclassical regimes for systems of quantum mechanics [Ll], [L7] and of quantum 



field theory |T2|, [Tj|, [19], |2(| |21j| . It makes use of the closed time path (CTP) functional 



technique, due to Schwinger and Keldysh [p2f. This is a path integral technique designed 



to obtain expectation values of field operators in a direct way [ 23[ . In the semiclassical 



regime, a tree level approximation is performed in the path integrals involving the system 
variables. In this approximation, the equation of motion for the expectation value of the 
system field operator is the semiclassical equation, which can be directly derived from the 
effective action of Feynman and Vernon [B, O, O, [TR ITS, EOfl. When computing this effective 



action perturbatively up to quadratic order in its variables, one usually finds some imaginary 
terms which do not contribute to the semiclassical equation. The key point of this approach 
is the formal identification of the contribution of such terms to the influence functional 
with the characteristic functional of a Gaussian stochastic source. Assuming that in the 
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semiclassical regime this stochastic source interacts with the system variables, and, thus, 
these become stochastic variables, equations of the Langevin type are derived for these 
variables. However, since this approach relies on a purely formal identification, doubts can 
be raised on the physical meaning of the derived equations. 

The second approach is based on the description of the transition from quantum to 
classical behavior in the framework of the consistent histories formulation of a quantum 
theory. The consistent histories formulation, proposed by Griffiths [23], and developed by 
Omnes |25j and by Gell-Mann and Hartle |26|, ^7|, was designed to deal with quantum 
closed (i.e., isolated) systems. It is thus believed to be an appropriate approach to quantum 
cosmology, where the quantum system is the whole universe. The main goal of this for- 
mulation is the study of the conditions under which a set of quantum mechanical variables 
become decoherent, which means that these variables can be described in a probabilistic 



way [^g, ^7|> |28|, [29], g(J. When the closed system consists on a distinguished subsystem 
(the "system", which is also often called an "open system") interacting with its environ- 
ment, Gell-Mann and Hartle proposed a mechanism for decoherence and classicalization 
of suitably coarse-grained system variables [26], 27]. This approach allows to evaluate the 
probability distribution functional associated to such decoherent variables and, under some 
approximations, to derive effective quasiclassical equations of motion of the Langevin type 
for such variables M, E7L E3, M B2 . 



In Sec. [TT] we show that that these two approaches can in fact be related. In this way, 
we see that, on the one hand, the second approach sheds light into the physical meaning of 
the first one. On the other hand, the first approach provides a tool for computing effective 
Langevin-type equations to the second one. A large portion of this section consists of re- 
formulating the mechanism for decoherence and classicalization of Gell-Mann and Hartle in 
the language of the CTP functional formalism. 

In Sec. |HT| , we use the results of this analysis to formally derive effective equations of 
motion for the gravitational field in a semiclassical regime. This derivation relies heavily on 
the results of the previous section. We find that, in the semiclassical regime, gravity might 
be described by a background metric, solution of the semiclassical Einstein equation, plus 
some stochastic metric perturbations. The equation for these perturbations, the semiclassical 
Einstein-Langevin equation, is seen to incorporate the effect of the lowest order matter stress- 
energy fluctuations on the gravitational field. 

In this paper we use the (+ + +) sign conventions and the abstract index notation of 



Ref. 33 , and we work in units in which c — Ti — 1. 
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II. Effective equations of motion from environment-in- 
duced classicalization 



A. The CTP functional formalism for a system-environment inter- 
action 



We start this section by sketching the CTP functional formalism applied to a system- 



environment interaction and its relation with the influence functional formalism of Feynman 
and Vernon M. For more detailed reviews of the CTP functional formalism, see Refs. |TB|. 



and for the influence functional formalism of Feynman and Vernon, see Refs. || |6|, [Fj, || ]||, 1G 



TT| , |T3| . For simplicity, we shall work in this section with a model of quantum mechanics, but 



all the formalism can also be formally applied to field theory. It is instructive to maintain 
in this section the explicit dependence on h. Let us consider a model of quantum mechanics 
which describe the interaction of two subsystems: one, called the "system", with coordinates 
q, and the other, called the "environment" , with coordinates We write the action for this 
model as S[q, Q] = S s [q) + S se [q, Q]%\ Let q{t) and Q(t) be the Heisenberg picture coordinate 
operators, which are assumed to be self-adjoint, i.e., q^ — q and = Q, and let q s and Q s 
be the corresponding Schrodinger picture operators. Suppose that we are only interested in 
describing the physical properties of system observables from some initial time ti until some 
final time tf > ti. Working in the Schrodinger picture, the state of the full system (i.e., 
system plus environment) at the initial time t = ti will be described by a density operator 
p s (ti). Let {\q,Q) s } be the basis of eigenstates of the operators q s and Q s . The matrix 
elements of the initial density operator in this basis will be written as p(q, Q; q' , Q'; ti) = 
s (q,Q\ p s {ti) \q',Q') s - For simplicity, we shall assume that the initial density operator can 
be factorized as p s (ti) =p s s {ti) <g> pl{U), in such a way that its matrix elements in coordinate 
representation can be written as p(q,Q;q',Q';ti) = p s (q,q';ti) p e (Q,Q';ti). However, the 
formalism can be generalized to the most general case of a non-factorizable initial density 
operator |36], [37|, p6| . We are interested in computing expectation values of operators related 
to the system variables only, for times t between ti and tf. The dynamics of the system in 
this sense can be completely characterized by the knowledge of the whole family of Green 
functions of the system. Working in the Heisenberg picture, these Green functions can be 
defined as expectation values of products of q(t) operators. These Green functions can be 
derived from a CTP generating functional in which only the system variables are coupled to 
external sources j+(t) and j-(t) || 0, [L3|, [H], [19|, |2D| . This CTP generating functional can 



^ven if, in order to simplify the notation, we do not write indices in these coordinates, q and Q have to 
be understood as representing an arbitrary number of degrees of freedom (which, in particular, can be an 
infinite number of degrees of freedom) . 

2 We shall assume that the action S[q, Q] is the one that appears in the path integral formulas for the 
model, which, in general, needs not to coincide with the classical action for the model p4l pq|. 
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be written as the following path integral^ 

ZV+J-] = fv[q + ] V[q._] p.{q +i ,q-.;t i )6{q +f -q- f ) ^ (sM^+nJdtj^-nfdtj^) ^ (21) 
with 

S c s[q+, q-\ = S s [q + ] - S s [q_] + S w [q+, (2.2) 

where Sip is the influence action of Feynman and Vernon, which is defined in terms of the 
influence functional JF IF as 

= jv[Q + ] V[Q_] Pe (Q +4 ,Q„ t ;t 4 ) 5(Q+ f -Q- f ) e i (^+,Q+]-^->Q-]). (2 . 3) 

We shall call S e s[q+,qJ\ the effective action of Feynman and Vernon. In these expressions 
we use the notation q+^ q+(ti), Q+f= <Z+(*/)> Q+i = Q+{ti)i Q+t = Q+(tf)-> an d similarly 
for q_ and Q_. All the integrals in t, including those that would define the actions S s [q] 
and S se [q, Q] in terms of the corresponding Lagrangians, have to be understood as integrals 
between t\ and tf. The CTP generating functional has the properties 

Z[j,j] = 1, Z[j_, U ] = Z*\j + ,jJ\, \Z\j + ,jJ\\< 1. (2.4) 

From this generating functional, we can derive the following Green functions for the system: 



T[q(t[)---q(t' s )]T[q(t 1 ).--q(t r )} 



iSj+ih) ■ ■ ■iSUitrK-i^j.it'J ■ ■ ■ (-i)5j_(t> s ) 



i 

i±=o 
(2.5) 

where h, . . . ,t r ,t[, . . . ,t' s are all between U and tf, T and T mean, respectively, time and 
anti-time ordering. The expectation value is taken in the Heisenberg picture state corre- 
sponding to the Schrodinger picture state described by p s (ti) at the initial time t = ti. The 
influence functional ( |2.3| ) can actually be interpreted as a CTP generating functional for 
quantum variables Q coupled to classical time-dependent sources q{t) through the action 
S se [q, Q] Let us consider the quantum theory for the variables Q in presence of classical 
sources q(t) corresponding to this action, and assume that the initial Schrodinger picture 
state for the quantum variables Q is described by the density operator p s e (U). For this the- 
ory, let U[q](t,t') be the unitary time-evolution operator, which can be formally written as 
U[q](t,t') — Texp -| //, dt"H s [q](t") , for t>f, where H s [q](t) is the Hamiltonian operator 



3 A way of generalizing the formalism to a non-factorizablc initial density operator consists in the fol- 
lowing |3(], ^6). One writes the initial density matrix in coordinate representation as p(q,Q;q',Q';ti) — 
p s {q,q'\ti) p se (q,Q;q' ,Q';ti), where p s is chosen in such a way that fdq p s (q,q;ti) = 1. Then, the CTP 
generating functional can be written as ( |2.l| ), with 

e | Sefi[q+,g- 



Jv[Q+]v[Q-] p i e(g+ i ,0+ i ;«- i ,Q- 1 ;ti)«(Q+ / -Q- / )e* (s[ * + ' Q+] - s[9 - g - ]) . 
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in the Schrodinger picture. This Hamiltonian operator depends on t as a function of q(t) 
and their derivatives q(t), and this gives a functional dependence on q in the operator IA. It 
is easy to see that [Z7|, 0, [TJ, 0, fffi 



^rF[g+,g-]=Tr ^(t i )Wt[g_](t / ,t i )W[ g+ ](t /) t i ) = (W* [<?_](*/, t f ) %+](*/, (2.6) 

where we use ( )^tt i ) to denote an expectation value in the state described by p s e {ti). From 
this expression, it follows that the influence functional satisfies 

•7tof[9>?] = ^iF[g-,g+] = ^ [<?+,<?-], \JTj?[q+,q-] \ < 1, (2.7) 

or, equivalently, in terms of the influence action, 

<SiF[g,g] = 0, Sw[q-,q+] = S£[q+,q-], lm.S w [q+, g_] > 0, (2.8) 

and similar properties follow for S , cff [g + , A decoherence functional for the system, where 
the environment variables have been completely integrated out, can now be introduced as 
the functional Fourier transform of the CTP generating functional in the external sources: 

Z\j+,j.] = Jv[q + ] 2%_] V[q + ,q_] e*/*^*-*-*-), (2.9) 
that is, from fl2.1|) we have that 

V[q + ,q_] = p s (q +v q_ v U) S(q +f - q. f ) e* s -*[*+'H (2 .1 ) 

In the consistent histories approach to quantum mechanics, V[q + ,q_] is known as the deco- 
herence functional for fine-grained histories of the system [|26], ^7], ^8], ^l|, |3^] . 

The environment of a system has to be understood as characterized by all the quantum 
degrees of freedom which can affect the dynamics of the system, but which are "not accessi- 
ble" in the observations of that system. This environment includes in general an "external" 
environment (variables representing other particles, or, in the context of field theory, other 
fields) and an "internal" environment (some degrees of freedom which, from the fundamental 
quantum theory point of view, would be associated to the same physical object as the "sys- 
tem" variables, but which are not directly probed in our observations of the system) 
For instance, a problem which has been studied using the influence functional method is 
that of quantum Brownian motion [|, g g 0, 0, 0, ||, 0, ||, 0, g|, HI HZL 0- In 



this problem, one is interested in the dynamics of a macroscopic particle interacting with a 
medium composed by a large number of other particles. In this example, one considers that 
the only "observable" system degree of freedom is the center of mass position of the macro- 
scopic particle, whereas the remaining microscopic degrees of freedom of the macroscopic 
particle are considered as environmental variables. Such "internal" environment degrees of 
freedom, and also those of the particles of the medium (the "external" environment), are 
usually modelized as an infinite set of harmonic oscillators. In the context of field theory, 
one would typically consider as "inaccessible" to the observations the modes of the field of 
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interest with characteristic momenta higher than some cut-off momentum JD], [T2|, [IB| . In 



the case of the gravitational field, this has been considered by Whelan [I2| in a toy model 
designed to investigate the decoherence mechanism for gravity. 

It is convenient at this stage to distinguish between these two kinds of environmental 
variables, so let Q represent the coordinates of the "external" environment (the coordinates 
of "other particles") and q v the "unobservable system" coordinates (the coordinates of the 
"internal" environment). As before, q will represent the "true" system coordinates. One 
could now simply replace Q by (Q, q v ) in the previous expressions. However, for conve- 
nience, we shall do the integrations in the environmental variables in two steps. The action 
of the full system will be now written as S[q, q v , Q], and, as before, we shall assume a to- 
tally factorizable initial density operator p s (t«) = p%ti) ® p^iU) <8> p s e {U)i which leads to 
an initial density matrix in coordinate representation of the form p(q,q v ,Q',q',q' v ,Q'',ti) = 
Psi.1i <?'; U) P u ( ( ?u' ?u> Pe(Q, Q'] U) (notice that we are now using the subindex e for the "ex- 
ternal" environment). Such a factorization is based on the assumption that the interactions 
between the three subsystems can be neglected for times t < tj. Unfortunately, in most 
situations, this assumption does not seem to be very physically reasonable, especially for the 
"true" system- "internal" environment interactions. One would need to consider the general- 
ization of the formalism to a non-factorizable initial density operator mentioned above and 
the analysis would be more complicated. We start defining 

e | (Sf[q+]-Sf [g_]+S s c ff[ g+ ,Q + ;g_,<2_] ) 

= /p[ 9u+ ]Pk_] Pu ( gu+i , gu ^ (2.11) 

where the effective action for the system S* [q] is chosen to be real and local. Notice that 
the effective action S^[q + , Q + ; g_, Q-] has analogous properties to those of Sw in (|2.8f) . We 
introduce now an effective influence functional and an effective influence action as 

J*§ [ q+ ,qJ\ = e m<^-] = Jv[Q + ] V[Q.] p^Q^Q^U) 6{Q +f - Q- f ) e i«[ g+> Q +;g -,Q-]. 

(2.12) 

With these definitions, the effective action of Feynman and Vernon, S e s[q+, q~], which ap- 
pears in expression (271) can be written as 

S eS [q+, g_] = Sf[q + ] - Sf[q.] + S$ [q + , q_]. (2.13) 

Note that, since S e s[q + ,qJ\ satisfies the same properties as S\p in ( p.8|) , it follows from the 
last expression that Sf§ has also these properties. 

B. The "naive" semiclassical approximation 

The usual "naive" semiclassical approximation for the system variables consists in performing 
a "tree level" approximation in the path integrals involving the q variables in expression ( |2.1| ) 
[|, 0, 0, [15|, [y| ^0| . Therefore, the CTP generating functional is approximated by 



J (S efl [qf [3},q { %]}+h$dtj + qf \j]-h f dt j 



Z\j + ,j-} * e B V— — - — - (2.14) 
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where q±[j] = q±[j + ,j-} are solutions of the classical equations of motion for the action 
S eS [q+,q-} +HJdtj + q + - hfdtj_q_, that is, 



5q ± {t) 



(2.15) 



which satisfy the boundary condition q+(tf) = g_ (i/)- Whenever this approximation is 
valid, we can see from (gjj) , (|2~T5|) and (U) that (g(t)) ~ q {0) (t), with g (0) = g™[j + =j_ = 
0] = q^ } [j+=j- = 0}, that is, g (0) (t) is a solution of the two equivalent equations: 



SS eS [q+,q-) 



5q + (t) 



0. 



5S eS [q+,q-] 



:,(«) 



-(*) 



0. 



(2.16) 



;,(») 



One can see that these two equations are actually the same equation, and that this equation 
is real. This is the semiclassical equation for the system variables. In a naive way, one would 
think that, when the above semiclassical approximation is valid, the system would behave as 
a classical system described by the coordinate functions q (0) (t), i.e., that one could substitute 
the description of the system in terms of the operators q(t) by a classical description in terms 
of the functions q (0) (t). However, one can see from (|2.14|) , (|2.15| ) and ( |2.5| ) that, in general, 

mt[) ■ ••<?«)] T[q(h) ■ ■■q(t r )\ > 4 Q (0) (h) ■ ■ ■ q m (t r ) q^(t[) ■ • V 0, (O- (2-17) 

Thus, in general, whenever the above approximations are valid, we can only interpret the 
solutions of the semiclassical equation as representing the expectation value of the operators 



C. Further coarse-graining and decoherence 

Decoherence takes place in a set of quantum-mechanical variables when the quantum inter- 
ference effects are (in general, approximately) suppressed in the description of the properties 
of a physical system which are associated to that variables. When this happens, such deco- 
herent variables can be described in an effective probabilistic way. In the Heisenberg picture, 
we will say that a set of variables decohere when the description in terms of the operators 
corresponding to these variables can be replaced by an effective description in terms of a 
set of classical random variables, in the sense that the quantum Green functions for such 
operators become approximately equal to the moments of the classical random variables. 
For the Green functions fl2.5|) , it is easy to see that this would hold in an exact way if the 
CTP generating functional (|2~T1) depended on the sources j± only as a functional $q[j + — jJ\ 
of the difference j + — j_, or, equivalently, if the decoherence functional ( |2.9|) could be writ- 
ten as V[q + ,q_] = V q [q + ] 5[q + — g_]. However, in practice, one finds that this condition is 



usually too strong to be satisfied, even in an approximate way |26|, ^7], |28|, |3T], [3^, |42 [. 
One needs to introduce further coarse-graining in the system degrees of freedom in order to 
achieve decoherence. Let us then introduce coarse-grained system operators, which corre- 
spond to imprecisely specified values of the system coordinates. In the Heisenberg picture, 
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such operators can be denned as 



g c (*) = £gP*(*)> 



(2.18) 



where Pq(t) is a set of projection operators, labeled by some variables q (these are often 
discrete variables) , of the form 

P q (t) = Jdqdq v dQ^(q-q) \q,q v ,Q,t)(q,q v ,Q,t\. (2.19) 

Here {\q, q v , Q, t)} is the basis of eigenstates of the operators q(t), q v {t) and Q(t), and 7 
is a real function. We shall assume coarse-grainings of characteristic sizes a, that is, such 
that the function ^(q — q) vanishes or has negligible values for q outside a cell Iq of sizes a 
centered around q. This means that 



j dq 7(9 - <?) f(q) - j i dq l{q - q) f(q), 



(2.20) 



for any function f(q). In addition, the function 7 must be chosen in such a way that the 
set of projection operators is (at least, approximately) exhaustive and mutually exclusive, 
which means that 



E^W 



P q (t)Pg>(t)=6 W P 9 (t), 



(2.21) 



where I is the identity operator. For specific examples of operators satisfying the above 



properties in an exact or in an approximate way, see Refs. |31 



Next, we can introduce a family of decoherence functions for coarse-grained histories of 
the system [^6], ^7|, ^8|. ^9|. ^0|. [JI], [32]. In order to do so, let us consider a set {t\, . . . , t^} of N 
instants of time, such that t/- < tk+i, k — 0, . . . , N, with to = U and t^+\ = tf. Introducing 
two sets of values of q associated to such set of instants, {<?+} = {q+ i: • • • , <1+ N } and = 
. . . , q~ N }, the decoherence function for this pair of "coarse-grained histories" of the 
system is defined as 



V c ({q+}, {?-}W... 



(ti,...,t N ) 



Tr 



P, 



(2.22) 



P q+1 (ti)pPq_ 1 (t 1 )---P q _ N (t N ) 

where p is the density operator describing the state of the entire system (system plus en- 
vironment) in the Heisenberg picture (T> c is often called decoherence "functional" in the 
literature, but, for each set {ti, . . . , tjv}, this is actually a function of 2N variables). These 
decoherence functions can be written in a path integral form as 

N 

Z>c({g+}, {?-})(*,...,*„) = V[q+]V[q.] ]lrtq + (t k )-q +k )rtq4t k )-q-. h )V[q + ,qJ\, (2.23) 

J k=i 

where T>[q + ,q_] is the decoherence functional for fine-grained histories of the system fl2.9Q . 
From the definition ( p.22j ) and the properties ( 2.21 ), one can show that these decoherence 
functions have the properties 

E E Z>c({?+}, {?-}) = 1, T> e ({q.}, {q + }) = Vt({q + }, {g_}), (2.24) 

{$+}{?-} 
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and that the diagonal elements of the decoherence functions (the values of those functions 
in the limit — > q +k ) are positive. For iV > 1, we can also see that, if we divide the 
set {ti, . . . , t N } into a subset of M < N instants, {t[, . . . , t' M } C {ti, . . . , t N }, with t[ < 
■ ■ < t' M , and the subset of the remaining L = M — N instants, denoted as {t", . . . ,t" L } [i.e., 
{t l5 ...,t N } = {t[, . . . , t' M } U {?(, . . . , t'i}}, then 

^({g+L^g-D^,...,^) = £ v c({q+} N ,{q-}N)(ti,-,tNh ( 2 - 25 ) 

{q+} L {v-} L 

with {g ± } M = {q±{t' x )i . . . , q±(t' M )}, {q±} L = {q±(t"), q±(f[)}, where we use the notation 
q±(tk)=q± k , for k=l, .. .,N, and {q ± } N = {q ±1 , ...,q ±N }. 

To make contact with the CTP formalism, let us introduce now, in analogy with fl2.9| ). 
a family of generating functions for the coarse-grained system degrees of freedom as the 
following Fourier series: 

Zc({j + }, {./ })./• t N ) = EE V c ({q + }, {q tjv) jEL&uH-u-i-nl-u), (2.26) 

{?+}{?-} 



where {j±} = {j± 1 , ■ ■ ■ ,j± N }- Note that the properties Q2.24Q for the decoherence functions 
are equivalent to 

Z c ({0}, {0}) = 1, Z e {{j.}, {j + }) = Z c *({j + }, (2.27) 

From the generating function fl2.26| ), we can compute the Green functions 



G7 mi n . ms (t'n ■ ■ -x r xi, ■ ■ -X) = (T[g c mi (^) ■ ■ -q? s X)\ m x X) ■ • -r «.)] 

with rii,..., n r , mi, . . . , m s G IN, {t[, . . . ,t' r } C {ti, . . . , t^} and {t'{, . . . , £"} C {ti, . . . , t N } 
(thus, r, s < N): 

r n r - n r 1.1 .1 .11 (-id) ni+ +n ' +mi+ +m " Z c ({j + }, {j-})(t u ...,t N ) 

cmi - mi l ' ■ ■ ■ ' r ' x ' ■ ■ ■ ' s) [dj+wr- ■ vhwri-dj-wr- ■ -h^-cop {J±}=W 

(2.28) 

where j±(tk) =j± k , for k — l,...,N. The property Q2.25 ) can also be written in terms of the 
corresponding generating functions as 

Z c {{j+} M , O'-Im)^,.,^) = z c ({j+} N , {i-L)(ti,..,^)| {i±}i={0} , (2.29) 

with the notation {j±} M = {i±(^) ; • • • ,J±(^m)}; an d similarly for {j±} L and {j±} N . Notice 
that this last property is consistent with (|2.28|) , in the sense that, for instance, G^ in2 (t' 1 ,t' 2 ) 
can be equally computed either from Z c ({j + } 2 , {j-}a)(^,f 2 )> or from Z c ({j+} N , {j-} N )(t!,...,t N ), 
with N >2. 

Having introduced the coarse-grained description of the system in terms of the operators 
q c {t), we can now sketch the decoherence mechanism for them. For the Green functions 
( |2.28| ), one can show that the decoherence condition described above holds in an exact way if 

9 



the generating function ( |2.26| ) depends on the sources j± k only as a function of the differences 
j+ k —j- k , i.e., as $q({j + — j-})(t lt ...,t N )- Then, introducing the Fourier series corresponding 



to <3>g, we can write 

Zc({j+},{j-})( tl ,...,t N ) 



®q({]+-3-})(t 1 ,-,t N ) = J2 V Q({<l})(tl,-,tN) e 



Y,k=l^{j+ k -3- k ) 



(2.30) 

Note from the last expression that, if we interpret the function Vq as the probability distri- 
bution for a set of random variables %, k—l,...,N, associated to the instants tk, then <I>g 
is the corresponding characteristic function. Therefore, from ( 2.28 ), we get 

•qn^iH Mv+miH \-m s ^ (Sa\\ 



c mi ■■■ m, 



5 t S J 



{i}={0} 



c' 



[dmr ■ ■ ■ [dj(f r )r [w?)p ■ • ■ pj(*z)i 

=EW)m») M • ■ -Pit) TO) • ■ -TO) = (9" l («i) • ■ ■9 nr (4) TO) ■ • -TO)), 
{<?} 

(2.31) 

where ( ) c means statistical average of the random variables, and we use the notation q(tk) = 
Qk, j(tk) = jk, for k = 1, . . . , N. Note that, if ( |2.30|) is satisfied, then the property (|2.29|) 
reduces to 

®q({j} N )(tl,.,t N ) 



'Pq{{^}u){A,-.,t' M ) 



{j} L ={0} ' 



(2.32) 
(2.33) 



or, equivalently, 

jvJ(ti,...,tjv)- 

This last property is a necessary condition for the probabilistic interpretation ( p. 31 ) to be 
consistent. 

The conditions for decoherence ( |2.30| ) can be written in terms of the corresponding de- 
coherence functions as 

N 

V c ({q + }, {q-}) (tu ...,t N ) = 'Pg{{q + }) {tl ,...,t N ) U 6 9+h9 _ h . (2.34) 

fc=i 

These are actually the conditions for decoherence of coarse-grained system variables as stated 
in the consistent histories formulation of quantum mechanics [BB, I 



28, pm KG, 31, 32 



Notice, from fl2.21|) , that ( 2.34j) is always satisfied for a single instant of time (i.e., when 

n=i) m. 

We can now check that the interpretation of Vq as a probability function is actually 
correct. From the second of the properties fl2.24| ), we have that V^({q}) = Vq({q}), i.e., 
Vq is real. Since the diagonal elements of the decoherence functions are positive, Vq({q}) 
is also positive. These two properties of 'P$({?})(t 1 ,...,t JV ), together with (|2.33|) , are enough 
to guarantee that it can be properly interpreted as the probability distribution for a set of 
random variables associated to the instants £i, . . . , t^. From the first of the relations ( |2.24|) , 
which yields J2{q} ^5 ({<?}) — 1 3 h follows that this probability distribution is normalized. 
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In practice, the conditions for decoherence described above will be usually only satisfied 
in an approximate way Approximate decoherence is typically achieved through a mechanism 
which was proposed by Gell-Mann and Hartle ETJ. To see how this works, note that, 
if we assume coarse-grainings of characteristic sizes a [see (|2.20|) 1, and using ( p.lOj ), we can 
write the decoherence function ( |2.23|) as 



V C ({Q + }, {*-})<*,.,*) - / 2>[gf ] V[q^] f[ 2>[g?>] V[q^} Ps (q^, g^; t t ) 6{q™ - q™) 

"W },{*»_} k=l 



N 



k=0 



(2.35) 



where each path integration / U[q±^], for k = 0, . . . , N, is over paths q±\t) with t G [ifc, 
being t = ti and tjv+i = tf, an d we have used a notation to indicate that these paths are 
restricted to pass through the cells Iq ± at the instants for k — 1, . . . , N. From (|2.13|) , the 
modulus of each factor exp(^ S e g[q+\ g_ ]) in the last expression is exp(— ilmS'fp [g+\ g^]). 
Then, if for every k = 0, . . . , N, Im Sfp[q+ , g_ ], which is always positive or zero, is much 
larger than ft whenever the differences | q^ — | are larger than some "cut-off" sizes 
the integrand in ( p. 35 ) will be only non-negligible for \q+^ — q^\ < S k \ If the charac- 
teristic sizes a of the coarse-graining satisfy a 3> S k \ then the "off-diagonal" elements of 
T> c ({q + }, {q-})(t ly ..,t N ) are negligible and one has approximate decoherence [f26], |2^j. We 
should stress that Sfp[g+,g_] is the result of integrating out both the "external" environ- 
ment degrees of freedom and also the system degrees of freedom which are "not accessible" 
to the observations (the "internal" environment). In general, these two integrations play an 
important role in the achievement of this sufficient condition for approximate decoherence. 
A characterization of the degree of approximate decoherence has been given in Ref. |3l| (see 
also Refs. |29|, |8|). 

Typically, d {k) can be estimated in terms of Atk = tk+\ — tk- When this is the case, one 
usually finds that the Gell-Mann and Hartle mechanism for approximate decoherence works 
provided all the time intervals satisfy At^ > At c , k = 0, . . . , N, where At c is sufficiently 
larger than some characteristic decoherence time scale t D (t D can be written in terms of a 
and some parameters characterizing the environment and the system-environment couplings) 
[ [2T1| , |27|, |Sl|. For At c one should take the smallest value compatible with a specified degree of 



approximate decoherence. In this sense, we can think of a coarse-graining as characterized 
both by the sizes a and by the time scale At c . 

D. Effective equations of motion for the system 

Assuming that the mechanism for approximate decoherence described in the previous sub- 
section works, an approximate effective description of the coarse-grained system variables 
in terms of a set of random variables [in the sense of Eq. ([2.31 )1 is available, at least for 
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instants of time satisfying At k > At c , for k = 0, . . . , N. The corresponding probability distri- 
bution Vq({q})( tl ,...,t N ) is given by the diagonal elements of the decoherence function ( |2.22p . 
We shall next make an estimation of this probability distribution. This follows essentially 
the derivation of Gell-Mann and Hartle in Refs. l26|, |27l|. For alternative derivations for 



more specific models, see Refs. |28|, ^2[ . Introducing the new variables q A = q+— g_ 
and g E = | (g + + g_), and similarly for q± k , and assuming that a ^> S k \ note first, from 



( |2.35|) , that the restrictions on the integration over q A coming from the coarse-graining can 



be neglected in the diagonal elements of this decoherence function. Therefore, using (|2.23 ) 
and Q2.10D , and writing S cS [q + , g_] = S cS [q A , g E ], we get 



P?({fe})(*i,..,t w ) - faq*) f[l 2 (q s (t k )-q* k )V f [qv], (2.36) 

where 



k=l 



V { [q s ] = J V[q A ] p s (g s . + §g A .,g E .-§g A .;t,) e^^l. (2.37) 
9 A (*/)=° 

At this stage, we introduce two simplifications in our analysis. First, we restrict our 
evaluation to coarse-grained system variables having significance only up to certain scales, 
larger enough than a, so that the random variables q k can be well approximated by con- 
tinuous random variables. This approximation can be implemented with the use of a set of 
approximate projection operators Pq(t), with q being continuous variables, which satisfy the 
properties (|2.21 ) in an approximate way (see Refs. 31, [32| for an example). Then, all the 



sums Y,{q} can be replaced by integrals JYlk=i d<jk and the functions Vq({q})(t 1} ,., t t N ) become 
probability densities. Second, as long as we are only interested in the dynamics of the system 
on time scales much larger than At c (At c is proportional to the decoherence time scale t D , 
which is typically extremely small, see Refs. p5| , |39| , |43| for some examples), we can take 
the continuous time limit in ( |2.36| ). In order to do so, consider the instants tk = U + kAt, 
k = 0, . . . , N+l, with At = (tf—ti)/ (N+l). Introducing functions q(t), such that q(tk) = <jk 
(assumed now to be continuous variables), and letting N— >oo in ( |2.36|) [replace g Efc by q~k], 
with (tf — ti) maintained finite (thus, At— >0), we get a probability distribution functional 
associated to some stochastic variables q{t) 0: 



/^fe]7 2 fe-9]^f[fe], (2-38) 



where j[q] is the functional corresponding to ElfcLi 7(<z(*fc)) m the limit ^^oo (some redef- 
initions in the parameters entering in the function 7(g) may be needed in order that such 



limit is well defined, see Refs. []32], ^] for an explicit example of how this limit is taken). 
Notice that, if we take the limit to the continuous in time and in the variables q k in ( |2.30|) , 
we get a functional $q[j] which is the functional Fourier transform of Vq[q}. Hence, <&q\j] can 
be interpreted as the characteristic functional for the stochastic variables q(t) 0. From the 
probability functional (2138) or, equivalently, from the associated characteristic functional 



[by functional derivation with respect to the sources j(t)], we can compute the Green func- 
tions G^ 1 r ^'^ m (t[, . . . , t' r ; t'(, . . . , t") with each of the instants in {t[, . . . ,t' r } being separated 
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from ti and from the remaining instants in this set by intervals larger enough than At c , and 
similarly for the instants in {t'(, . . . , t"}. 

We can get a good approximation to the path integral ( f2.37[ ) by expanding [q A , q s ] 
in powers of q A and neglecting higher than quadratic terms, i.e., we make a Gaussian ap- 
proximation in this path integral. This expansion can be made using ( |2.13| ) and writing 



S , fp[g + ,g_] = S'fpfg^gJ. In this expansion, the dependence of S c fi[q A ,q s ] on the velocities 
g A (t)[] (we assume that there is no dependence on time derivatives of higher order) gives rise, 
after integration by parts, to boundary terms proportional to q A{ (we use that q A ^ = 0). For 
instance, assuming that 5'| ff [g] = Jdt L s (q(t),q(t),t), in the expansion of the terms Sf 1 we 
find a boundary term —Psiq-z^qx^t^q^, where p s = dL s /dq are the canonical momenta. 
Similarly, if Sfp depends on q A (t), its expansion will contain some boundary terms. However, 
since, in general, Sfp depends non-locally on q A (t) and q^(t), these terms will be more com- 
plicated. Note that we are considering models slightly more general than the ones studied 



by Gell-Mann and Hartle in Refs. p6|, |27 |, since we allow for the possibility of an influence 
action depending on q A (t) and q s (t). The motivation for considering such a generalization 
is that we are interested in field theory actions with interaction terms depending on the 
derivatives of the fields. 

One can show that, when expanding up to quadratic order in q A , the general form for the 
boundary terms in S§ is -Fi[g s ](£j) q A + iF 2 [g E ](tj) q\. + ifdtq A (t)F 3 [q s ](t,ti)q A , where 
Fi, F 2 and F 3 are real functionals of g E , which vanish when Sfp does not depend on q A (t). 
Finally, we get the following expansion: 



S cS [q A , g E ] = Sf + | q A ] - Sf [g E -i q A ] + S$[q A , q s ] = -ft [&](*.) q Ai + iF^U) q\. 

+ ijdt q A (t) F 3 [g E ] (t, U) q A . +Jdt q A (t)C[q s } (t) + ^-J dt dt' q A (t) q A (t') C 2 [&](t, t') + C>(g A ) , 

(2.39) 

with 

x qcS r i 

Pi[q)(U) = pMuM + Fi[q](U), C[q]{t) = + C x [q\(t), (2.40) 



and 



we*, t^i^ 1 sts ^ 



Sq A (h) ■ ■ ■ 5q A (t k 



(2.41) 



9 A =° 



where the functional derivatives with respect to q(t) are defined for variations which keep 
the value of q(t) fixed at t — ti and t = tf. 

Substituting the expansion (|2~39| ) into Eq. ( ET37D , we get a Gaussian path integral which 



can be calculated. Note that, since ImSff > 0, C 2 [q){t,f) is positive semi-definite. In 
order that the Gaussian approximation that we have carried out is valid, we must assume in 
addition that C 2 [g](i, t') is strictly positive definite and, thus, detC 2 [g] ^ 0. We get 

P t [q] ~ NWi[q] [det(C 2 [g]/27r^)] _1/2 e-^ d ' d ' ,C7 MW c "M(''*') c M(*'), (2.42) 



We understand that a term depends on q&(t) if it does so before any integration by parts. 
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where N is a normalization constant, C 2 1 is the inverse of C 2 defined by 

dt"C 2 (t,t")C 2 1 (t",t') = 5(t-t'), 



(2.43) 



Wi[q] = W{q{ti),p[q]{ti),U[q]{U)]U), with 

dq 



and 



W(q,p,U;U) e-^e-ho n p s ( qi + l %Ai _ 1^), ( 2 .44) 



p[q] (U) = Pl [q] (U) + h dt dt' F 3 [q] (t,U) C, 1 [q] (t, t') C[q] (f) 



n[g](ti) = F 2 [g](U) - -jdtdt'FM{t,U)C 2 x [q]{t,-t')FM{t'A)- (2.45) 



Note that the function W defined in (|2.44|) is a generalization of the Wigner function asso- 
ciated to the initial state of the system, and it reduces to the ordinary Wigner function for 
II = [4"4"|. Note that, in expression ( |2.42| ), the momenta p[q]{ti) in this generalized Wigner 
function are in general different from the canonical momenta p s (qi,qi,ti)- I n the case of Sfp 
non-depending on the velocities g A (t), one has p[q]{ti) = p s (qi,qi,U) and n[g](tj) = 0, thus, 
Wi[q] is the standard Wigner function. From the definition (|2.41|) , and using the properties 
of Sfp[q+, qJ\, we can see that 



Ci [?](*) 
Ca[?](t,0 



S Re Sip [q+, q 



5q+(t) 



^Sf§[q + ,q- 



h 

2 



-(*) 



q+=q-=q 

5 2 Im Sf§ [q+ , q_] 5 2 Im Sf§ [q+ , q. 



5q + {t)5q + {t') 
and then, from (|2.40|) and (|2.13 ), we have 



6q + (t)5q4t') 



<?+=<?- 



i+=i- 



(2.46) 



C[q](t) 



5S eS [q+,q-] 



Sq+(t) 



(2.47) 



Substituting Q2.42 ) into ( 2.38|) , we see that the only non-negligible contribution to the 
path integral in ( [2.38P come from those paths which are not very far deviated from the paths 
q {0) (t) which satisfy C[q <0) }(t) = 0, that is, which satisfy the semiclassical equation ( |2.16| ). 
This implies that only those paths q(t) which remain always near from the semiclassical paths 
q (0) (t) will give a non-negligible value to Vq[q}. In this sense, the mechanism proposed by 
Gell-Mann and Hartle is a mechanism for decoherence and classicalization of coarse-grained 
system variables. However, we see that, in general, Vq[q] has a complicated functional 
dependence on q(t). 

Let us then study the deviations from a specific solution of the semiclassical equation, 
that is, we shall now restrict our considerations to those paths q(t) which are distributed 
around a given solution q (0) (t) of the semiclassical equation. We can now introduce stochastic 
variables Aq(t) = q(t) — q m (t) which describe the deviations from q {0) (t). The associated 
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probability distribution functional V Aq [Aq] is equal to Vq[q (0) + Aq] up to a normalization 
factor, which, from ( J2.38 ) , is given by 



V q [q {0) + Aq] ~ J V[q] 7 2 [g] Vf [q m +Aq + q). (2.48) 

In practice, it is difficult to work out the explicit dependence of the probability distribu- 
tion functional on the characteristic parameters of the coarse-graining, a and At c , even in 
simple models [ ^BL Nevertheless, if such parameters are small enough so that the values 
of Vf [q m + Aq + q] do not change very much for the different paths q(t) which give a non- 
negligible contribution in (|2.48|) , the functional (|2.48j ) can be approximated by Vf [q {0) + Aq]. 



We can make a further approximation by expanding Pf[g (0) + Aq] around g (0) . This can be 
done by setting q s = q (0) + Aq in (|2.39|) , expanding in Aq, and substituting the result for 



this expansion in ( p.37|) . The result to lowest non-trivial order is 



V Aq [Aq] ~ N[q m ] Wi[q m + Aq] e-^^^'^^+AglWC,- 1 ^]^')^^)^]^^ ( 2 .49) 

where iV[g (0) ] is a normalization factor and Cx[g (0) + Aq] is the expansion of C[q {0) + Aq] 
to linear order in Aq. Notice that, in this probability functional, the factor Wi [g (0> + Aq] 
contains all the contribution arising from the initial state of the system. This generalized 
Wigner function, even if computed expanding around g (0) , will have in general a complicated 
non-local dependence on Aq, except when Sfjf is independent of q A , in which case it reduces 
to the standard Wigner function for the initial state of the system and depends only on Aqi 
and Aqi. If the deviations from g (0) are small enough, we can approximate Wi[g Co) + Aq] ~ 
iyj[g (0) ]. Then, with these approximations, the variables Aq are distributed in such a way 
that C i,[q {0) + Aq]{t) are Gaussian stochastic variables characterized by 

(C L [q m + Aq](t)) c = 0, (C L [q^+Aq](t)C L [q^+Aq](t')) c = C 2 [q^]{t,t>). (2.50) 

Thus, the equation of motion for Aq is the Langevin equation 

C L [ q ™+Aq](t)+m = 0, (2-51) 
where £(t) is a Gaussian stochastic source with 

<£(f)> c = 0, mm) c = C 2 [q m ](t,t'). (2.52) 

We should mention that there are very simple models for quantum Brownian motion in 
which all the actions involved are quadratic in their variables and the interaction terms are 
independent of the velocities ||, [5|, |^, [U], ^6|, [27], [28], [31], |37|, |40| . For such models, assuming 



that the environment is in an initial state of thermal equilibrium, the influence functional 
can be computed exactly and it is Gaussian. The effective action of Feynman and Vernon 
in these cases is exactly of the form Q2.39p , with Cf[q s ](t) linear in g E , C 2 (t,t') independent 



of q s and Fi = F 2 = F 3 = 0. Thus, for these models, expression ( |2.42| ) is actually exact. 
In these cases, with the approximation Vg[q] ~ Vf[q], one can derive a Langevin equation 
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for the stochastic variables q(t), without need of introducing a specific solution g (0) of the 
semiclassical equation. This Langevin equation is simply C[q](t)+£(t) = 0, being £(£) a 
Gaussian stochastic source with (£(£)) c = and (£(t) £(t')) c = C 2 (t,t')- However, for models 
with more complicated actions, we are only able to derive effective equations of motion for 
the deviations Aq around a given solution q {0) of the semiclassical equation. 



E. A quick method to obtain the Langevin equation 

Starting with the effective action of Feynman and Vernon (|2.13| ), there is a quick way to 
obtain the Langevin equation (|2.51| ) for the deviations Aq around a specific solution of the 
semiclassical equation. This method has actually been extensively used in the literature, 
in the context of quantum Brownian motion |], [II], [17J, and also in the context of field 
theory |I^, |I^, [TJ], ^0], [21], including some models for gravity interacting with a scalar field 
[@j B 0' HH> EH' EH • One starts with an expansion of this effective action around a solution 
q (0) (t) of the semiclassical equation up to quadratic order in perturbations Aq± satisfying 
Aq + (ti) = Ag_(tj) and Aq + (tf) = Ag_(t/) (in the simplest models, in which this effective 
action is exactly quadratic in q + and g_, one works directly with the exact expression). From 
( |2.39|) , it is easy to see that the expansion for the influence action reads 

S$ir + Aq + ,r + Aq_}= |^(Ag + (t)-Ag_(t))C 1 [g (0) + i(Ag + + Ag_)](t) 

+ k I dtdt ' ( A 1+(t)-^Q-(t)) C 2 [q m ](t,t') (Ag + (i')-Ag-(*')) + 0(Ag 3 ), (2.53) 

where it is understood that C\ has to be expanded up to linear order. Using the identity, 
which follows from a Gaussian path integration, 

g-^j / dtdt' (Aq + (t)-Aq-(t)) C 2 [g(°)](t,t') (A q+ (t')-Aq-(t')) _ J p^j e iJdtS(t) (Aq+(t)-Aq-(t)) ^ 

(2.54) 

where V^] is the Gaussian probability distribution functional for the Gaussian stochastic 
variables £(t) characterized by ( |2.52| ), that is, 

e -i/didi'€(t)c 2 -v o) ] (t,t')m 

we can write in this approximation 

[g (0) +Ag+,g (0) +Ag_]| = e -k^Sf§[ q ^+A q+>q ^+A q .) = / e ±fdH(t)(A q+ (t)-Aq-(t))\ ? ( 2 .56) 



where ( ) c means statistical average over the stochastic variables £(t). Thus, the effect of 
the imaginary part of the influence action ( |2.53| ) on the corresponding influence functional is 
equivalent to the averaged effect of the stochastic source £ (t) coupled linearly to the pertur- 
bations Aq± (note that, in the above expressions, the perturbations Ag ± are deterministic 
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functions). Notice that expression ( |2.54| ) or, equivalently, ( |2.56| ) give the characteristic func- 
tional of the stochastic variables [[J. The influence functional, in the approximation 
( J2.53 ) , can then be written as an statistical average over £: 

F§ [g (0) +Ag + ,g (0) +Ag_] = (e* ^f^^l) , (2.57) 



where 

^f[Ag + ,Ag_;£] = Re ^ e F ff [g (0) + Ag + , g (0) + Ag_] +/dt£(f) (Ag+(i) - Ag_(i)) + 0(Ag 3 ), 

(2.58) 

where Re S'fp can be read from expression fl2.53|) . The Langevin equation ( |2.51|) can be easily 
derived from the action 

A ff [Ag + , Ag_; £] = [g (0) +Ag + ] - Sf [g (0) + Ag_] + A$[Aq+, A?_; £], (2.59) 

where 5| [g Co) +Ag-|_] has to be expanded up to second order in the perturbations Aq±. That 



is. 



<L4eff[Ag + ,Ag_;£] 



<5Ag + (i) 



(2.60) 



A(j + =Ai}_=Ai3 



leads to Eq. ( g3g ). 



III. Effective equations of motion for the gravitational 
field 

In this section, we shall apply the results of the previous section to derive effective equations 
of motion for the gravitational field in a semiclassical regime. In order to do so, we will 
consider the simplest case of a linear real scalar field $ coupled to the gravitational field. 
We shall restrict ourselves to the case of fields defined on a globally hyperbolic manifold A4. 
In this case, we would consider the metric field g a b{x) as the system degrees of freedom, and 
the scalar field and also some "high- momentum" gravitational modes, considered as 
inaccessible to the observations, as the environment variables. Unfortunately, since the form 
of a complete quantum theory of gravity interacting with matter is unknown, we do not 
know what these "high-momentum" gravitational modes are. Such a fundamental quantum 
theory might not even be a field theory, in which case the metric and scalar fields would 
not be fundamental objects. Thus, in this case, we cannot attempt to evaluate the effective 
actions in Eq. ( |2.11| ) starting from the fundamental quantum theory and integrating out 
the "high-momentum" gravitational modes . What we can do instead is to adopt the usual 
procedure when dealing with an effective quantum field theory. That is, we shall take for the 
actions Sf ff [<7] and S^[g + , $ + ; g~, $_] the most general local form compatible with general 
covariance and with the properties of [these properties are analogous to those of Sif in 
Eq. (|2T8|)1 |35| , f45f . The general form for S'| ff [^] is 

1 (R - 2K B ) + a B C abcd C abcd + (3 B R 2 + ■■], (3.1) 



sf[g] 



d x 



-g 



16nG B 
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where R and C a bcd are, respectively, the scalar curvature and the Weyl tensor associated 
to the metric g a b, 1/Gb, A b /Gb, O-b and [3b are bare coupling constants and the dots 
represent terms of higher order in the curvature [because of the Gauss-Bonnet theorem in four 
spacetime dimensions, there is no need of considering terms of second order in the curvature 
different from those written in Eq. (ft.l|) 1- Since M. is a globally hyperbolic manifold, we can 
foliate it by a family of Cauchy hypersurfaces E t , labeled by a time coordinate t. We use 
the notation x for spatial coordinates on each of these hypersurfaces, and ti and tf for some 
initial and final times, respectively The integration domain for all the action terms must 
now be understood as a compact region U of the manifold Ai, bounded by the hypersurfaces 
S 4i and E t/ (i.e., as in the previous section, integrals in t are integrals between ti and tf). 
For the matter part of the effective action, let us consider the following ansatz: 



where £ is a dimensionless coupling parameter of the field to the scalar curvature, and the dots 
stand for terms of higher order in the curvature and in the number of derivatives of the scalar 
field. Self-interaction terms for the scalar field could also be included but, for simplicity, we 
shall ignore them in this paper. One can see that general covariance and the properties of 
5|g [g + , $ + ; g~~, $_] imply that imaginary terms and terms mixing the "plus" and "minus" 
fields in this action must be necessarily non-local. Thus, within a local approximation, 
the ansatz ( |3.2[ ) is the most general form for this action. We shall comment below some 
limitations of this local approximation. 

In order to simplify the analysis, we neglect the contributions of the higher order terms 
not written in Eqs. ( |3.1| ) and (|3.3|) . Assuming that the mass of the scalar field is much smaller 
than the Planck mass, this is a good approximation in a regime where all the characteristic 
curvature scales are far enough from the Planck scales. The terms in the gravitational 
Lagrangian density proportional to R 2 and C a b c dC abcd need to be considered in order to 
renormalize the matter one-loop ultraviolet divergencies. 

Assuming the form ( |3.2| ) for the matter part of the effective action, we can now introduce 
the corresponding effective influence functional as in Eq. ( |2.12| ). Let us assume that the 
state of the scalar field in the Schrodinger picture at the initial time t = ti is described by 
a density operator p s (tj) (in the notation of the previous section, this was pKU), but, here, 
we drop the subindex e to simplify the notation). If we now consider the theory of a scalar 
field quantized in a classical background spacetime (A4, g a b) through the action (|3.3| ), to this 
state it would correspond a state in the Heisenberg picture described by a density operator 
p[g}. Let {|<p(x)) s } be the basis of eigenstates of the Schrodinger picture scalar field operator 
<f> s (x): $ s (x) \ip) s = <^(x) \tp) s . The matrix elements of p s (ti) in this basis will be written as 
Pi[ip, 0] = s ((p\ p s {ti) \<p) s - We can now introduce the effective influence functional as 



(3.4) 



[g + , $ + ; g-, $_] = S m [g + , $ + ] - S m [g- , $_] 



(3.2) 



with 




(3.3) 
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and the effective influence action will be given by J r fp[g + ,g~] = e lS ^\s + >9 \ 

Of course, trying to show how the mechanism for decoherence and classicalization of 
the previous section can work in this case would involve some technical difficulties, such as 
introducing diffeomorphism invariant coarse-grainings and eliminating properly the gauge 
redundancy (with the use of some suitable Faddeev-Popov method) in the path integrals. 
We are not going to deal with such issues in this paper. We shall rather assume that they 
can be suitably implemented without changing the main results for the effective equations 
of motion. 



Expression ( |3.4| ) is actually formal, it is ill-defined and must be regularized in order to 
get a meaningful quantity for the influence functional. We shall formally assume that we 
can regularize it using dimensional regularization, that is, that we can give sense to Eq. Q3.4| ) 
by dimensional continuation of all the quantities that appear in this expression. We should 
mention that, however, when performing specific calculations, the dimensional regularization 
procedure may not be in all cases the most suitable one. In this sense, one should under- 
stand the following derivation as being formal. Using dimensional regularization, we must 
substitute the action S m in (|3.4|) by some generalization to n spacetime dimensions. This 
can be taken as 

S m [g, ®r] = -\j d n x ^ [g ab d a <5> n d b <S> n + (m 2 + $*] , (3.5) 

where we use a notation in which we write a subindex n in all the quantities that have 
different physical dimensions than the corresponding physical quantities in the spacetime 
of four dimensions. The quantities that do not carry a subindex n have the same physical 
dimensions than the corresponding ones in four spacetime dimensions, although they should 
not be confused with such physical quantities. A quantity with a subindex n can always be 
associated to another one without a subindex n; these are related by some mass scale /i, for 
instance, it is easy to see that <J> n = [i 2 <p. 

In order to write the effective equations for the metric field in dimensional regulariza- 
tion, we need to substitute the action (|3.1| ) by some suitable generalization to n spacetime 



dimensions. We take 

1 2 



Sf[g}=^- 4 Jd n xV^9 



(R — 2Ab) + - OLE [RabcdR a ° — RabR a ) + @bR 



1QhG b v ' 3 

(3.6) 

where R a bcd is the Riemann tensor and, again, the mass parameter /i has been introduced 
in order to get the correct physical dimensions. Using the Gauss-Bonnet theorem in four 
spacetime dimensions, one can see that the action obtained by setting n = 4 in ( |3.6| ) is 
equivalent to ( |3.1| ). The form of the action ( |3~^| ) is suggested from the Schwinger-DeWitt 
analysis of the divergencies in the stress-energy tensor in dimensional regularization [46 |. 
The effective action of Feynman and Vernon (|2.13|) is in our case given by S c s[g + , g~] = 



^s S [d + ] ~ Ss S [d~] + Sjf [s f+ ' 9~\- Since the action terms (|3.5| ) and ( |3.6|) contain second order 
derivatives of the metric, one should also add some boundary terms to them [?]]. The 
effect of these boundary terms is simply to cancel out the boundary terms that appear when 
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taking variations of S e g[g + , g~] that keep the value of g£ b and g~ b fixed on the boundary of U. 
They guarantee that we can obtain an expansion for S c a[g + , g~] analogous to ( p2.39| ) , with no 
extra boundary terms coming from the integration by parts of terms containing second order 
derivatives of g£ b = g£ b — g^- Alternatively, in order to obtain the effective equations for the 
metric [equations analogous to ( |2.16|) and (|2.51 )1, we can work with the action terms ( |3.5|) 
and ( p.6|) (without boundary terms) and neglect all boundary terms when taking variations 
with respect to g^ b . From now on, all the functional derivatives with respect to the metric 
must be understood in this sense. 



A. The semiclassical Einstein equation 

From the action ( |3.5| ), we can define the stress-energy tensor functional in the usual way 

T ab ia, *„](*) = 6S r [9 :V , (3.7) 

yj-g{x) 5g ab {x) 

which yields 

r'b^J = v «$ n V^n-^ n V c *«Vc^-^ a6 m 2 ^+( - V"V k + G ni ) $^ (3.8) 

where xj a is the covariant derivative associated to the metric g ab , □ = VaV a > an d G ab is 
the Einstein tensor. Working in the Heisenberg picture, we can now formally introduce the 
stress-energy tensor operator for a scalar field quantized in a classical spacetime background, 
regularized using dimensional regularization, as 

f: b [g] = T ab [g, f ab [g] = /i-(™- 4 ) f^[g], (3.9) 

where & n [g] (x) is the Heisenberg picture field operator in n spacetime dimensions, which 
satisfies the Klein-Gordon equation 

(□ - m 2 - £R) <f> n = 0, (3.10) 

and where we use a symmetrical ordering (Weyl ordering) prescription for the operators. 
Using Eq. ( |3.10| ), one can write the stress-energy operator in the following way: 

f: b [g] = ~{s7 a $n[g], V h ®n[g]} + V ah [g]& n [gl (3.11) 

where T> ab [g] is the differential operator 

Vf =(c~i) 9 ab {x)U x + £ (R ab (x) - vJvi) , (3-12) 
being R ab the Ricci tensor. From the definitions (|3.4|), ( |3.7| ) and Q3.9|), one can see that 

2 8S§[g + ,g- 



-9\ x ) 



$9ab( x ) 



(T?(x))[g], (3.13) 



g^=g~=9 
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where the expectation value is taken in the n-dimensional spacetime generalization of the 
state described by p[g]. 

As in Eq. ( [2. lb] ), if we derive S cS [g + , g~] with respect to g+ b and then set g+ b = g~ b = g ab , 
we get the semiclassical Einstein equation in dimensional regularization: 

G ab [g] + A B g ab ) - {% a B D ab + 2(3 B B ab ) [g] = (T n a %], (3.14) 



8%G B 

where the tensors D ab and B ab are defined as 
1 8 



D ab _ 



ab 



J d n x V^g (Rc def R cdef - R cd R cd ) = \ g ab (R cdef R cdef - R cd R cd + or) 

- 2R acde R b cde - 2R acbd R cd + AR ac R b - 3 UR ab + V a \J b R, (3-15) 



and 

B ab = -j=jjj- J dn x V^gR 2 = g 9 abR2 ~ 2RRab + 2 V a V b R - 2g ab nR. (3.16) 

From equation ( |3.14j ), after renormalizing the coupling constants in order to eliminate the 
divergencies in /i~^ n ~ 4 ^ (T£ b ) [g] in the limit n^A and then taking this limit, we will get the 
semiclassical Einstein equation in the physical spacetime of four dimensions: 

G ab [g] + Ag ab ) - 2 (aA ab + (3B ab ) [g] = (f a R b )[g]. (3.17) 



8ttG 



In the last equation l/G, A/G, a and (3 are renormalized coupling constants, (T^lg] is 
the renormalized expectation value of the stress-energy tensor operator, and we have used 
that, for n = 4, D ab = (3/2) A ab , being A ab the local curvature tensor obtained by functional 
derivation with respect to the metric of the action term corresponding to the Lagrangian 
density C abcd C abcd . 



B. The semiclassical Einstein-Langevin equation 

According to the results of the previous section, assuming that some suitably coarse-grained 
metric field satisfies the conditions for approximate decoherence and that the approximations 
of subsection [II][D] are valid in a certain regime, small deviations from a given solution g ab of 
the semiclassical Einstein equation (|3.17|) can be described by linear stochastic perturbations 



h a h to that semiclassical metric. These perturbations satisfy a Langevin equation of the form 



( 2.51|) , which shall be called the semiclassical Einstein-Langevin equation. Our next step will 
be to write the semiclassical Einstein-Langevin equation in dimensional regularization. Let 
us assume that g a b is a solution of Eq. ( |3.14j ) in n spacetime dimensions. The semiclassical 
Einstein-Langevin equation in dimensional regularization has then the form 

Gf[g+h] + A B (g ab -h ab )) - ( | a B D? + 2{3 B B L b ) [g + h] = (f?\ [g + h] 



8ttG B \ L 1 V ' / \3 L ^ B i J IV' I ■ \J- n j L \ 

+ 2fi~^C b , (3-18) 
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where h ab is a linear stochastic perturbation to g a b, h ab =g ac g bd h c d, that is, g ab — h ab + 0(h 2 ) 
is the inverse of the metric g a b + h a b, and, as in the previous section, we use a subindex l to 
denote an expansion up to linear order in h ab . In this equation, (T^ b )[g + h] is the expectation 
value of T^ b [g + h] in the n-dimensional spacetime generalization of the state described by 
p[g + h], and £™ is a Gaussian stochastic tensor characterized by the correlators 



0. 



abed \ 



with [see Eqs. ( p32|) and (glgfl ] 



2 A^p] (*,</) 



g(x)J-g(y) 



6 2 lmS^[g + ,g-} S 2 lmSf§[g + ,g 



(3.19) 



(3.20) 



9^=9 =9 
a6\ 



We can write Eq. ( |3.i£| ) in a more explicit way by working out the expansion (T% °) [g-\-h\. 
Since, from Eq. ( |3.13 ), we have that 



8S§[g+h + ,g+h- 



[g + h] 



det(g + h)(x) 



$htb(x) 



(3.21) 



h+=h~ 



this expansion can be obtained from an expansion of the influence action Sfp [g + h + ,g + h~] 
up to second order in h^ b (in this expansion, we can neglect boundary terms). At the same 
time, we can obtain a more explicit expression for the noise kernel ( j3.20| ). To perform this 
expansion for the influence action, we have to compute the first and second order functional 
derivatives of Sjp [g + , g~] and then set g^ b = g~ b = g ab . If we do so using the path integral 
representation ( |3.4| ), we can interpret these derivatives as expectation values of operators 
in the Heisenberg picture for a scalar field quantized in a classical spacetime background 
(A4, g ab ) as, for instance, in expression (|3.13|) . The relevant second order derivatives are 



S 2 S§[g\g- 



y/-g(x)y/-g(y) 6 gtb{ x ¥ gtd{y) 
1 S 2 S§[g + ,g-) 



-H: bcd [g](x,y)-KT d [g}(x,y) + l N 1 



abed r 



jabed 



y/-g(x)y/-g(y) 5 glb( x )Sg C d(y) 

with 



9"*~=9~ =9 



g^=g~=g 



HZ bcd {g}(x,y)- l N« bcd {g}(x,y), 



abed r 



x,y), 
(3.22) 



N: 



abed 



^y) = i(\ T n\*) - (T:\x)) , T: d (y) - (T n c %)) \)[g], 



-tab 



~-\cd I 



*CM*,v) = ^{T\T: b (x)T: d (y)))[g], 



H: bcd [g](x,y) = - 



- /- 

4 \2 



K 



abed 



x,y) 



f: b (x),f^(y)}J[g], 
1 / 5 2 S m [g^ n ] 



yj-g(x)yj-g(y) \Sgab(x)5g cd (y) 



(3.23) 
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using again a symmetrical ordering (Weyl ordering) prescription for the operators in the last 
of these expressions. All the expectation values in these expressions are in the n-dimensional 
spacetime generalization of the state described by p[g]. In the above equations, { , } and 
[ , ] mean, respectively, the anticommutator and the commutator, and we use the symbol 
T* to denote that, first, we have to time order the field operators and then apply the 
derivative operators that appear in each term of the product T ab (x)T cd (y), where T ab is the 
functional (13. 81). For instance, 



T*( V ;^)v^n(^ 

(3.24) 

where T is the usual time ordering. Notice that all the kernels that appear in expressions 
( |3~22D are real. 

In fact, from ( |3.23|) , we see that the noise kernel N^ bcd , and also the kernel H£ bcd , are 
free of ultraviolet divergencies in the limit n — > 4. This is because, for a linear quan- 
tum field, the ultraviolet divergencies in (T^ b (x)T^ d (y) X J are the same ones as those of 

T£ b (x)) (T^ d (y)Y Therefore, in the semiclassical Einstein-Langevin equation 



can perform exactly the same renormalization procedure as the one for the semiclassical 
Einstein equation (|3.14j ). After this renormalization procedure, Eq. ( |3.18|) will yield the 
semiclassical Einstein-Langevin equation in the physical spacetime (n = 4). It can be written 

as 

_L( G ^ + h ]+A (g ab - h ab )) - 2 (aAf + 0B*) [g + h] = (f*) L \g + h] + (3.25) 
being £ afe is a Gaussian stochastic tensor with 

(C\x)) c = 0, (C\x)C d (y)) c = N abcd [g}(x,y), (3.26) 

where N abcd = lim n ^ 4 ^~ 2{n ~ A) Nf cd . Notice from ( ^23l) that the noise kernel N abcd [g](x, y) 
gives a measure of the lowest order fluctuations of the scalar field stress-energy tensor around 
its expectation value. Thus, the stochastic metric perturbations h ab , solution of the semi- 
classical Einstein-Langevin equation ( p.25| ), account for the back reaction of such matter 
stress-energy fluctuations on the spacetime geometry. For a more detailed analysis of the 
semiclassical Einstein-Langevin equation and some of its applications, see Ref. P7[ . 

Going back to the expressions in dimensional regularization, which may be useful for 
calculational purposes, we can now write the expansion of the influence action around a 
given metric g ab . From ( |3.13|) and (|3.22|) , taking into account that Sfp[g,g] = and that 
s§ [g-,g + ] = [g + , g~], we get 



Sg\g+h + ,g+h-] = -Jd n x^g~[x){f:\x))[g} (K b {x)-K b {, 
~J d n xd n y sf^f^y) (H^[g](x, y)+K a n b " 1 [g](x, y)) (h+(x)h+(y)-K b (x)h; d (y) 
-\ fd n xd n y^^)sPgV)Hf n cd [g}(x,y) (K b (x)h^(y)-h- b (x)ht d (y) 
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+ l - [d n xd n yJ^)f^)K bcd [g](x,y) (h$(x)-h*(x)) (K d (y) ~ KM) + 0(h 3 



(3.27) 



From ( p.23[ ), it is easy to see that the kernels satisfy the symmetry relations 



H« bcd (x, y) = H s cdab (y, x), H« bcd (x, y) = ~H cdab {y, x), K a n bcd (x } y) = K c n dab (y, x). 



Sn 



Using these relations, and defining 



H a n bcd (x, y) = H^ d (x, y) + H^ d (x, y), 



(3.28) 
(3.29) 



we can write the expansion (|3.27|) as 



S^[g+h+,g+h-) = -Jd n x^-g(x) {TT{x))[g} [h ab {x)\ 

-\fd n x d n y f^ix) yf^jj) M*)l ( H n cd [9) y) +K bcd [g] y)) {KM} 

+ l - J d n x d n y yf^)yf^) [h ab (x)\ K bcd [g](x, y) [h cd (y)] + O(^), (3.30) 
where we have used the notation 

[hab} = K b -h~ b , {h ab } = h+ b +h~ b . (3.31) 
Using this expansion and noting, from (|3.23|) , that 

1 1 / ST ab \n <b .1/Vl \ 

M 



K**\g](x y) = - 1 - (f*(x))\g]^^S»(x- V ) - 1 _J_/ *^ **](*) 

(3.32) 

we get, from ( |3.21f) , 

(t\x)) L [g+h] = (t\x))[g] + (f^ b [g; h](x))[g] - 2 J d n y f^y) H^[g](x,y) h cd (y), 

(3.33) 

where the operator T^ 1)ab is defined from the term of first order in the expansion T£ [g+h, $ n ] 



as 



T»\g+h, *„] = T ab [ 9l $ n ] + T^ ab [g, $ n ; h], T^ ab [g; h] ee T^ ab [g, $ n [g\; % (3.34) 

using, as always, a Weyl ordering prescription for the operators in the last definition. Note 
that the third term in the right hand side of Eq. (|3.33|) is due to the dependence on h cd of the 
field operator $ n [g + h] and of the dimensional regularized version of the density operator 
p[g + h}. 



Substituting ( |3.33| ) into (|3.18j) , and taking into account that g ab satisfies the semiclassical 
Einstein equation ( 3.14Q , we can write the Einstein-Langevin equation (|3.18| ) as 



(G^ ab [g; h](x)-A B h ab (x)) - \ a B D^ ab [g; h](x) - 2(3 B B^ ab [g; h](x) 



^(f^ ab [g; h}(x))[g] + 2 d n y f^y) ^^H^[g\{x yV ) h cd (y) = 2^\^{x). 



(3.35) 
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In the last equation we have used the superindex (1) to denote the terms of first order in the 
expansions Gf[g + h], Df\g + h\ and Bf[g + h\. Thus, for instance, Gf[g + h] = G ab [g] + 
G (1)ab [g;h}. The explicit expressions for the tensors T Wab [g, $ n ; h], G Wab [g;h], D {1)ab [g;h] 
and B (1)ab [g; h] are given in the Appendix. From T {1)ab [g,$ n ; h), we can write an explicit 
expression for the operator T^ 1)ab . Using the Klein-Gordon equation (|3.10|) , and expressions 
(|3. 11 ) and ( 3.12|) for the stress-energy operator, we can write this operator as 

Ti 1)ab [9; h] = C- g ab h cd - 6 a c h b d - 6X) f^[g] + T ah [g- h] & n [g], (3.36) 



where T \g\ h) is the differential operator 



ab 



- g ab h c : 

2 c 



-I 



V c V a ^ + V c V b K - Oh ab - \7 a V b K - g ab X7 c V d h cd 



ab , 



+ g ab oh c c + (y a h b c + v b K - Vch ab - 2g ab V d h cd + g ab \7ch d d ) XJ C -g ab h cd V c V c 



(3.37) 



and it is understood that indices are raised with the background inverse metric g ab and 
that all the covariant derivatives are associated to the metric g ao . Substituting expression 
(|3.36|) into Eq. ( |3.35| ), and using the semiclassical equation ( |3.14| ) to get an expression for 



y. 4 ) (T£ b ) [g] , we can finally write the semiclassical Einstein-Langevin equation in dimen- 
sional regularization as 



1 



G 



(l)ab 



8nG B 
1 
2 

_ -(n-A) jrab(Q2 



-9 ab G cd h cd 



G ac h b + G bc K + A R (h 



-,ab 



■-g ab h' 
2 ' 



\u B [D^ b 



-- g ab D cd h cd + D ac h b c + D bc h a c )(x) - 2p B (B Wab -^g 



1 
2 

abed 



ab B cd h cd + B ac h b c + B bc h a c ) (x) 



[g] + 2j d n y f^j^ n ^HZ bcd [g] (x, y) h cd (y) 



2fi 



-{n— 4) tab 



(3.38) 



where the tensors G ab , D ab and B ab are computed from the semiclassical metric g a b, and 
where we have omitted the functional dependence on g a b and h a b in G Wab , D Wab , B (1)ab and 
T ab to simplify the notation. Notice that, in Eq. Q3.38|) , all the ultraviolet divergencies in 
the limit n— >4, which shall be removed by renormalization of the coupling constants, are in 



$^(x)y and the symmetric part H^ b (x,y) of the kernel H° (x,y), whereas, as we have 
pointed out above, the kernels N® (x, y) and H£ bcd (x, y) are free of ultraviolet divergencies. 
Once we have performed such a renormalization procedure, setting n = 4 in this equation 
will yield the physical semiclassical Einstein-Langevin equation, Eq. (|3.25| ). Note that, due 
to the presence of the kernel H% bcd (x,y) in Eq. (|3.38|) , such Einstein-Langevin equation will 
be non-local in the metric perturbation. 



C. Discussion 

We have seen that effective equations of motion for the metric field of the form ( |3.17| ) 
and ( |3.25| ) follow from the local approximation ( |3.2| ) for the effective action describing the 
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"effective interaction" of the metric and the scalar field. A more realistic evaluation of 
this effective action starting from a fundamental theory of quantum gravity would certainly 
lead to some real and imaginary non-local terms in this action. In some situations, the 
contribution of these terms to the effective equations of motion for the metric (note that they 
would also give some extra terms in the semiclassical equation) might not be negligible and, 
in any case, one would expect that their role in the decoherence mechanism for the metric 
field would be important. This would represent non trivial effects coming from the "high- 
momentum" modes of quantum gravity, which are not part of the gravitational field described 
by the classical stochastic metric g ab + h ab , but which can be source of this gravitational 
field in the same way as the matter fields. The contribution of these neglected terms to 
the equations for the background metric g ab and for the stochastic metric perturbation h ab 
would be similar to the contribution of the scalar field through its stress-energy operator, but 
with this operator replaced with some "effective" stress-energy operator of such primordial 
"high-momentum" gravitational modes coupled to the scalar field. These equations would 
take the form (|3.17 ) and ( |3.25|) only when the effect of this "effective" stress-energy tensor 
on the classical spacetime geometry can be neglected. A way of partially modelizing this 
effect would consist on replacing the stress-energy operator Tf[g] by Tf[g] + tf[g], where 
is the stress-energy tensor of gravitons quantized in classical spacetime background 
(M,g ab ) 0. 

We end this paper with some comments on the relation between the semiclassical Einstein- 
Langevin equation Q3.25|) and the Langevin-type equations for stochastic metric perturba- 
tions recently derived in the literature || |7|, ||, [L3|, [U], [15], [Tj|. In these previous derivations, 



one starts with the influence functional (|3.4j ), with the state of the scalar field assumed to 
be an "in" vacuum or an "in" thermal state, and computes explicitly the expansion for the 
corresponding influence action around a specific metric background. One then applies the 
method of subsection [II][E| to derive a Langevin equation for the perturbations to this back- 
ground. As we have seen in subsection [TI|[E[ this method yields the same equations as the one 
used in this section. However, in most of the previous derivations, one starts with a "mini- 
superspace" model and, thus, the metric perturbations are assumed from the beginning to 
have a restrictive form. In those cases, the derived Langevin equations do not correspond 
exactly to our equation, Eq. (|3.25| ), but to a "reduced" version of this equation, in which only 
some components of the noise kernel in Eq. ( |3.26|) (or some particular combinations of them) 
influence the dynamics of the metric perturbations. Only those equations which have been 
derived starting from a completely general form for the metric perturbations are actually 
particular cases, computed explicitly, of the semiclassical Einstein-Langevin equation (|3.25| ) 

m Bill- 
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Appendix: Expansions around a background metric 

For a metric of the form g a b = g a b + h a b, where h a b is a small perturbation to a background 
metric g ab , we list the expansions of metric functionals around the background metric up to 
linear order in the perturbation. In the following expressions, all the tilded quantities refer 
to functionals constructed with the metric g a b, whereas that the analogous untilded ones are 
constructed with the background metric g a b- In particular, y a and \/ a are respectively the 
covariant derivatives associated to the metric g a b and to the metric g a b, and y a = g ab \/b) 
□ = y a Va> V a = 9 ab \/b, 1=1 = V a Va> where g ab and g ab are respectively the inverses of g a b 
and g a b- We shall also raise indices in the metric perturbation with the inverse background 
metric g ab : h a h = g ac h cb and h ab = g ac g bd h cd . 

g ab = g ab -h ab + 0(h 2 ), (A.l) 

~g = ^g~(i + \K + o(h 2 )y (A.2) 
K b = T a b + \ (VaK + v b K - v c M + o(h 2 ). (A.3) 

For a scalar function /, 

VaVbf = VaVbf ~ \ V7 (VAc + VbKc ~ VAb) + 0{h 2 ), (A.4) 

□/=□/- V a vV K b - vY (V b h ab - \ Vahl) + 0(h 2 ), (A.5) 

v a vY = vV/ - v a v c f h b c - vV7 K-\ v7 (v a h b c + v b K - v c h ab ) + o(h 2 ). (A.6) 

For a tensor t ab , 

at ab = at ab - Xj c Xj d t ab h cd + {g ae \jH db + g be \jH ad - ± g cd xj e t ab ) {xj c h de + V Ae - VeM 

+ ^ (<T^ + 9 bc t ad ) (v e Vd h ce + D/iod - V e Vc ^e) + 0(/* 2 ), (A.7) 
For the curvature tensors, 

Rab = Rab + ^ (V C Vahc + V°Vb Kc ~ &Kb ~ Va Vb K) + 0(> 2 ), (A.8) 

R a b =R a b - R c b K + l - (v c VbK + v c V a h c - ah a b - VbV a K) + 0(h 2 ), (A.9) 

R = R- R ab h ab + V a V b h ab - Oh a a + 0(h 2 ), (A.10) 
l ab = R ab -R ac h b c -R bc h a c - 
G ab = G ab + G Wab + 0(h 2 ), with 

QWab = _ R ac h b_ R bc h a + 1 [ R h ab + g ab R cd fl cd +V C V^c + VV^c " D/l" 6 ~ V° V^c 

+ g ab (nh c c - vVM] = -G ac h b c - G bc h a c + X - [-Rh ab + g ab R cd h cd 
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R ab = R ab -R ac h b c -R bc h a c + \ (v°V a h b c +V c V b K - Oh ab - V a V b K) + (^ 2 ), ( A - n ) 



c 



+V c V a ^c +V c V b K - ah ab -\7 a V b K + g ab (ah c c - vVM], (A.12) 

R\cd = R a bcd + \ (VcVbK + VcVdK + VdV a ^ c - VcV a h d - VdVbK 

-V d VcK) + 0(h 2 ), (A.13) 

Rabcd = Rabcd + - {R e bcd^ae + Ra cd^be) + ~ (VcVb^ad + VdVa^&c ~ VcVa^&d 

- VrfV^a C ) + 0(/i 2 ), (A. 14) 

j^abcd j^abcd ^ ^c^j^abce j^d _|_ ^j^abed j^c _|_ j^aecd j^b _|_ j^ebcdj^a^ 

+ \ {v c V b h ad + \7 d V a h bc - \7 c V a h bd - y d y b h ac ) + 0(/i 2 ), (A.15) 

v a v 6 # = v a v^ - vVOR cd M + v fl vV c v d ^ - v a v b °K - v a v c R K 

- V V^R K~\ V C R (V a h b c + X7 b K - X7ch ab ) + 0(h 2 ), (A.16) 
DR = DR-n (R ab h ab ) + □ V a V'hab ~ 2 h a a - V a V b R h ab 

- V a R (V b hab - \ Vah b b ) + 0(h 2 ), (A.17) 

nR ab = UR ab - a(R ac h b c + R bc h a c ) - S7 c V d R ab h cd - \7 c R ab (V d h cd - \ Vch d d ) 
+ V c R ad (Vch b d + Vdh b c - v'M + x7 c R bd (VcK + sj d h a c - y a h cd ) 

+ \ R ac (v d Vch b d + uh b c - v VM + ± R bc (v d VcK + uh a c - v VM 

+ 1 (□ y c V a ^c + D V c V b K - a 2 h ab - □ y a V b K) + 0(h 2 ), (A.18) 

R 2 = R 2 - 2RR ab h ab + 2Ry a xj b h ab - 2R Uh a a + 0(/i 2 ), (A.19) 
^^ ab = RR ab - i?i? ac ^ - RR bc h a c - R ab R cd h cd + ^R {xj c \j a h b c + \j c \j b h a c 

- nh ab - S7 a V b K) + R ab (vV^W - Oh c c ) + 0(h 2 ), (A.20) 
R ab R ab = R ab R ab - 2R ab R c a h bc + R ab (2 V c vAc - Uh ab - VaVbK) + 0(h 2 ), (A.21) 

R ac R\ = R ac R b - R ac R bd h cd - R cd (R a c h b d + R b c h a d ) + ± R ac (xj d XJ c h b d + SJ d SJ b h cd 

-Oh b -VcV b h d d ) + \ R bc (V d VcK + V d V a h cd - Uh a c - VcV «h d d ) + 0(h 2 ), (A.22) 

R abCd Rabc d = R^^Rabccl ~ 2R abcd ' R a bce h d + ^R ahcd \J c \J b h ad + 0(ll 2 ), (A. 23) 

^>acM^> ^ R acbd R d -\- — R d (^R acde fob _|_ J^bcde y^a c^j^acbe yd c^j^bcae fod^ 

2 

+ \ R acM (S7 e Vch de + V e Vdhce - Ohcd - VcVdK) - \ R cd (2 sj c y d h ab 

+ V a V^ cd + V b V a h cd - 2 y a V c ^ d - 2 V b V c h ad ) + 0(/i 2 ), (A.24) 

ryacde r)b jyacde -nb ^ / r>acde r>f ib , rjbcde r>f ua \ r> -nacde r>b if 

K K cde~ K K cde ~ 2 ^ K cde> l f + K K cde h f ) ~ ZK K cdf h e 

+ \ R^ iVdVcht + Ve - Ve " VrfV^ce) 
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+ ^ R bcde (VdVcK + VeV^cd - VeVcK - X7dV a h ce ) + 0(h 2 ), 



(A.25) 



B ab = B ab + B mab + 0(h 2 ), with 

5 (1)afe = -i i? 2 /i ab - g ab RR cd h cd + 2R {R ac h b c + R bc h a c ) + Rxj a xj b h c c + 2R ab {R cd h cd + Uh c c 

- V VM + g ab V c V d R h cd + 2URh ab -2 \j c \j a R h b c - 2 V c V b R h a c + g ab V V(^M 

- v c [R (g ab Vch d d + v a h b c + v b K - Vch ab )} + 2 9 ab a (R cd h cd + uh c c - v VM 

-2 V a V b (R cd h cd + Uh c c - vVM, (A.26) 
D ab = D ab + D (1)ab + 0{h 2 ), with 

D (i)ab = 1 (tfrffl^ _ R cdef Rcdef } h ab + (flacdefcV + ^cde^a/ ) _ ^ ^jjacfcM 

+ 4 J R 6c /i ad + i? acde /i^ + i? 6cde /i^ - 2i? acbe ^ - 2R bcae h d e ) + # ab (i? c/ ,R cp - R cdef R cdeg ) h 9 f 

- AR ac R bd h cd + 4R acde R b cdf h{ + ^R cd {2 y c y d h ab + V a V^«j + vV°^ 



- 2 V a Vc^S " 2 V b Vc^S) + 2 ^ (vV^S - 7D^ + 7 ^ " 4 VcV & ^) 

+ ^ R^ (V d VcK ~70h a c +7 V^hcd ~ 4 VcV"^) " ^ (VdVc/^e + VeV^cd 

- VdV^ce - VeVc^d) - (VdVc^e + VeV^cd ~ VdV^ce ~ VeVcK) 

- l - g ab R cd (2 V e Vc^e - - VcV^:) + 2^^ Ve V<A/ - R acM (V e Vch de 
+ V e V^ ce - ah cd - sjcSjaK) + i V c ^ (Vch ab - v a h b c - v b K) - 3 v c R ad (Vch b d 



Vdh b c - v b Kd) - 3 v c R bd (VcK + VdK- v a M - \ (g ab v c R - 6 v c R ab ) (2 v d K, 



- V ch d d ) ~\°R h ab - X7 a V c R h b c - X7 b V°R K~\ 9 ab V°V d R hcd + 3 vV^ h cd 

+ 1 □ {2R ac h b c + 2R bc h a c + uh ab + v a V b K - V c V a h b c - V c V b K) — ^ g ab n (i^/icd + D/i 

- vVM - vV(# c %d + Oh c c - S7 c S7 d h, 



cd , 



(A.27) 



For the stress-energy tensor functional, 

T ab [g, $„] = V a $n V 6 $n - ^ <? ab V c $n Vc$n - ^ ^m 2 $ 2 + f (^D - V " v' +G a6 ) $ 2 , 
T a ^, $ n ] = T tt6 [^, $ n ] + T« a6 [(7, $ n ; fe] + 0(/i 2 ), with 



T<^ fe [^, $ n ; h\ = -T ac [g, $„] ^ - T bc [(?, $„] ^ - - ( V c iVci + ™H 2 ) ^ 6 

+ i ^ a6 v c $ n V d $n ^cd + I [— -R h ab + g ab R cd h cd + v°V a h b c + V c V b K - V° V 6 ^ 
- n/^ + ^ (n^ - vVM + (V a ^c + y b K - S7ch ab - 2(? a6 S7 d hcd + ^ ab Vc/id) V c 



+ 2h ab U-2g ab h cd x 7 c x 7 



ab ; 



c 



$ 2 . 



(A.28) 
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